Quantum back-action in spinor condensate magnetometry 
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We provide a theoretical treatment of the quantum backaction of Larmor frequency measurements 
on a spinor Bose-Einstein condensate by an off-resonant light field. Two main results are presented; 
the first is a "quantum jump" operator description that reflects the abrupt change in the spin state 
of the atoms when a single photon is counted at a photodiode. The second is the derivation of a 
conditional stochastic master equation relating the evolution of the condensate density matrix to the 
measurement record. We comment on applications of this formalism to metrology and many-body 
studies. 
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Atomic vapor magnetometers of spin-polarized alkali 
atoms are among the most sensitive field sensors demon- 
strated to date p. These magnetometers, based on 
the optical detection of Larmor precession, have demon- 
strated field sensitivities in the attoTesla/Hz 1 / 2 regime 
[3J- The use of optically trapped ultracold atoms as the 
sensing medium holds promise for magnetic microscopy 
at high spatial resolution as well as for significant im- 
provements in field sensitivity via entanglement-assisted 
techniques [3H5]. Spinor Bose-Einstein condensates are 
particularly suited to field sensing applications due to 
their low spin relaxation rates and absence of density- 
dependent collision shifts [5J. 

The detection of Larmor precession and the subsequent 
estimate of the magnetic field relies on the dispersive in- 
teraction between the collective atomic spin and the op- 
tical field, followed by a quantum-limited measurement 
of the light. This interaction entangles the optical and 
atomic degrees of freedom. Measuring the light breaks 
this entanglement and must therefore cause a backac- 
tion on the atomic spins. The ultimate sensitivity of 
such atomic magnetometers is governed by the interplay 
between the projection noise of the atomic spins, pho- 
ton shot noise and the quantum backaction due to the 
measurement of the optical field. Here, we provide a 
theoretical treatment of this backaction by means of a 
conditional stochastic master equation that relates the 
condensate evolution to the optical measurement record. 
One of the main results that distinguishes this work from 
similar studies, e.g. |7], is the derivation of a quantum 
jump operator reflecting the backaction on the atoms in- 
duced by the detection of a single photon. 

Model: We consider a magnetometer consisting of a 
spinor Bose-Einstein condensate of n 3> 1 spin-1 bosons 
trapped in their motional ground state. Our formalism is 
well suited to both spinor gases with ferromagnetic (e.g. 
87 Rb) and those with polar (e.g. 23 Na) ground states. 
We assume in this paper that the spatial extent of the 
gas is small enough that the single-mode approximation 



(SMA) is valid. We also note that our results can be 
readily adapted to systems with different internal state 
spaces or expanded to incorporate spatial variations in 
the atomic spin field. 

The atoms are prepared initially with high fidelity in 
the spin state, 

H) =c + | + l)+c |0}+ C _|-l) (1) 

in e.g. the z basis, with the normalization ^ i \ci\ 2 = 1. 

Before proceeding further, we introduce a "vector" no- 
tation over the spin indices for clarity and compactness, 
with the definitions 

v = (v + ,v ,v-), (2) 

|v| = v + +vo+V-, (3) 

v! = v + \v lv-\, (4) 

v w = (5) 

and the dot product defined in the usual way. For later 
use, we also define probabilities 

Pi=\<k\ 2 , (6) 

where the state normalization implies that |p| = 1. With 
this notation, the full n-particle state can be compactly 
expressed as 

|* n ) = aJ l |o)/VSi= Y, \ t c ^ h (7) 

IJI=n ' J - 

where aK — ca} is the Schrodinger field creation operator 
for state |0) is the vacuum, and |j) is the state with 
ji atoms in spin eigenstate i. 

The spin state of the condensate is optically detected 
either by measuring the phase imprinted onto a short, 
coherent pulse of the light field by the condensate as in 
phase-contrast imaging 8, 9 , or by monitoring the polar- 
ization of the light field as in polarization spectroscopy 
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(see, for example [TO])- Without loss of generality, we 
assume the former method via balanced homodyne de- 
tection. 

By introducing a general vector of coupling strengths 
g, we can describe the light-atom interaction via the 
Hamiitonian 



The signal in Eq. ( 11 ) can be interpreted as follows. If 



Hi = hg ■ Nb r b. 



(8) 



where A, — a\di and b is the light field annihilation oper- 
ator. To derive this Hamiitonian, we have assumed that 
the probe field is detuned sufficiently far from the atomic 
resonance so that the excited atomic states can be adi- 
abatically eliminated [TT]. The effect of this Hamiitonian 
is to imprint a phase rotation on the optical field that 
contains information about the spin state of the atoms. 
The optical transit time is taken to be short enough that 
other contributions to the condensate's evolution may be 
neglected for now. The quantum state of the combined 
atomic spin and two output optical fields after passing 
through the final beam splitter, just before measurement, 
is therefore 



^(l^ e - lG 'J))®|^(l- le -* G -j)), (9) 



where Aq is the dimensionless amplitude of the coherent 
optical pulse (i.e., the mean photon number is I^IqI 2 )? t is 
the light-matter interaction time and G = rg. The probe 
beam's contribution to the output fields has been phase 
rotated through the unitary evolution given by the Ham- 
iitonian ([8]). Importantly, unless the atoms are in the un- 
likely configuration of a Fock state of the spin field, this 
interaction entangles the optical field with the collective 
spin state of the atoms [12] . In the rest of this letter, we 
quantify the backaction induced on the condensate due 
to measurement of the light. 

Assuming the photodectors have unit efficiency, the 
joint probability of obtaining photon counts C+ and C_ 
at the first and second photodiodes, respectively, is 




C+!C_! 



C++C- 



Ul=n 



x (l + sinG-j) c +(l-sinG-j) c '-.(10) 



In the case n\G\ <C 1 we find for the mean observed 
photon count difference the simple expression 



(AC) 



|A | 2 nG-p 

\A a \ 2 n (g + G+ ~ G - (F, 



G+ — 2G + G- 2 , 
2 {F*. 



(11) 



where AC = C+ — C_ is the difference between photon 
clicks at the first and second photodiodes, and (F z ) = 
p + — p- and (F^) — p + + p_ are the single atom spin 
expectation values. 



the atoms are precessing in a magnetic field perpendicu- 
lar to the light, the phase contrast signal has a large DC 
component proportional to the number of atoms, an AC 
signal at the Larmor precession frequency, and a third, 
smaller signal at twice the Larmor frequency. The am- 
plitude of the respective contributions can be controlled 
by changing the polarization and/or the detuning of the 
probe light, thereby changing the relative values of the 
quantities G [8]. 

In an experiment, the number of condensed atoms is 
classically uncertain, so the spin state is replaced 

by a density matrix p. Atom number in an experimental 
condensate obeys roughly Poissonian statistics: 



TJ 



(12) 



n=0 



where n is the mean atom number over many runs and 
\^ n ) is given by Eq. Q. Equation (111 is modified for 



this case by making the substitution n — > n. 

So far, we have considered the measurement of the con- 
densate state by a pulse of light. To reach a continuous 
measurement limit, we instead consider coherent light 
with a photon flux rate per unit time / rather than an 
amplitude Aq, and examine the evolution of the system 
in a time interval St such that f6t < 1. In this case, 
the probability of counting no photons at either detec- 
tor is nearly unity, and any case where multiple photons 
are counted is negligible. By making the substitution 
Aq — > v / fSt we can therefore model this continuous ob- 
servation as a series of short, weak pulses. The free evo- 
lution of the atoms will be negligible on this timescale, 
i.e., between individual photon clicks. 

In this limit, the detection of a photon by either detec- 
tor (and no photon by the other) has the effect of acting 
on the density matrix with one of two jump operators: 



P(1,0) 
/3(0,1) 



J+PiJ+, 



J+ = 



(13) 



The traces of these post-click density matrices represent 
the probabilities of the respective photodetections. Both 
probabilities are much less than one; under our assump- 
tions, no click (and hence, no change except a small 
rescaling of p) is by far the most likely outcome. In the 
standard theory of photodetection applied to a cavity 
with damping k and single mode operator a, the jump 
operator for a photon click is \J nSta 13 . In the present 
case, there are two major differences: First, no atoms are 
removed from the system, and hence the jump commutes 
with N. Second, as long as (G-N) <C 1, a single jump will 
have only a very small effect on the atomic state. This 
is to be expected - a single non-resonant photon pass- 
ing through our apparatus ought not perturb the atoms 
much. In the next section, we combine this jump opera- 
tor with the free spin evolution of the atoms to arrive at 
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a master equation describing the effects of the detection 
process. 

Conditional stochastic master equation: In accordance 
with the experimental situation, we consider both the 
atom-light interaction and the Zeeman interaction in- 
duced by the ambient magnetic field, 



H z = %lMbB • F, 



(14) 



where <7l is the Lande g-factor, /ib is the Bohr magneton, 
B is the applied magnetic field - eventually the field to 
be detected - and F is the Schrodinger field spin opera- 
tor. This Hamiltonian could be modified or replaced in 
more general settings, e.g. by inclusion of the quadratic 
Zeeman effect. We consider the dynamics of the system 
over a timescale At 3> St and an incident flux / such 
that a natural separation of scales occurs, 



fAt 

g L p B (B-F)At 
(G-N) 



(15) 
(16) 
(17) 



Crudely speaking, these scales posit a large photon num- 
ber, small free evolution, and an even smaller effect from 
a single quantum jump, respectively. 

Next, we consider the set of all possible photocounts 
{(G+, G_)} and the probabilities of their detection over 
the interval At. The density matrix will be acted on 
by jumps J± alternated with unitary evolution due to 
Eq. ( 14 ) . Because we do not have access to the micro- 



scopic details of photon arrival order nor photon arrival 
times, we must sum (integrate) over the final density 
matrices for all possible arrival orders (times) consistent 
with a particular photocount for the interval At. This 
calculation is aided by the fact that the jump opera- 
tors will commute with the unitary evolution operators at 
least to order e 2 . The resulting conditional density matri- 
ces are normalized by the probability of detection of their 
particular photocounts. See [11] for further details or [14] 
for a thorough treatment of a similar situation. We use 
the resulting probability distribution on {(C+, C_)}, and 
find that, to next to leading order in e, 



v+ 



fAt 



(1±(G-N)), 



(18) 



that is, the mean equals the variance in the counts. The 
covariance between G+ and G_ vanishes to this order. 
Also, we note that the means reproduce the signal seen 
in Eq. (111. Because the higher moments of the distri- 



bution agree with the equivalent Gaussian moments to 
next to leading order in e, we are able to characterize the 
photocounts as independent Gaussian random variables, 



G± = |(l±(G-N»At+y|(l±(G-N»AW±, (19) 
where AW± are independent Wiener increments, with 
(AW|)b = At, 



Here the subscript E refers to the ensemble average (that 
is, over many experimental runs) rather than the usual 
expectation value. These Wiener increments are Gaus- 
sian random variables of zero mean that are functions of 
time. In the continuum limit, A — > d, these equalities 
are exact without ensemble averaging. 

The result Eq. ( 19 1 is crucial to the derivation of the 



complete stochastic master equation. We apply the oper- 
ators J± each C± times to p, act on it by the free unitary 
evolution, divide by its trace, and take the continuum 
limit to obtain |11) . 

f = \\P, H 1 ] + ^ (G • Np + pG • N - 2(G • N)p)£_ (t), 

(21) 

where the modified Hamiltonian H' is 



H' = h[g L p B B-F+-G-N[f 



/£+(*) 



(22) 



dV± 



are uncorrelated Gaussian white noise 



and £±(t) - dt 
functions related to the photocurrents, with dV± 
(dW+±dW-)/V2. 



Equations ( |21[ ) and ( 22 ) succinctly encapsulate the full 
effects of the measurement scheme on the atoms. First, 
the Hamiltonian is modified by a steady, effective inter- 
action resulting from Eq. ds|, with the photon number 
operator replaced by a classical coherent flux. Shot noise 
in the photon flux leads to an additional random, but 
unitary, evolution. Finally, the second term in (21) is 



the non-unitary evolution induced by the measurement 
of G • N. 

Examples: We first examine a very simple application 
of this master equation to illustrate its salient features 
in an analytically closed form before moving on to the 
more experimentally relevant case of the measurement 
of a magnetic field. Namely, we will assume that there 
is no external field, B = 0, and coupling coefficients 
G+ = — G_ = G, Go = 0. This means that the setup 
will measure only F z , i.e., 



G • N = GF Z 



(23) 



and the total Hamiltonian is quantum non-demolition 
(QND) for the observable F z . 

We can derive equations of motion for the observables 
of the system by multiplying Eq. (21) by their opera- 
tors and taking the trace. Because of the large num- 
ber of atoms involved, we assume that the Gaussian 
state ansatz will be valid; thus, only equations for the 
observables, their variances, and their covariances are 
needed. The chain rule of Ito calculus (resulting from 
Eq. (20|) must be used when considering the deriva- 



tives of the (co)variances, as they are non-linear func- 
tions of the other variables. We obtain, for example, 
d(F z )/dt — GyJ 'f£-(t)v z , where v z is the variance of (F z ), 
and v z = —fG 2 v z , which quickly integrates to 



(AW+AW-) E = 0. 



(20) 



v z (t) = K(o)- 1 + fGhy 



(24) 
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FIG. 1: The evolution of the condensate's mean spin along 
each axis for I = 0.017 sa t, normalized by atom number, for a 
single simulated experimental run. The optical detection of 
Larmor precession has virtually no effect on the free dynamics. 
Insert: The normalized photocurrent difference oscillates at 
the correct frequency. 



Note that determination of the mean value depends on 
the experimental record through £-(t), whereas the vari- 
ance always decreases as the measurement continues. 
This should be unsurprising - continuous observation of 
F z leads to greater certainty in its value. Over long times, 
the uncertainty can drop below the standard quantum 
limit. Next, we introduce a magnetic field to be mea- 
sured, and present the results of a numerical integration 
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FIG. 2: The evolution of the condensate's mean spin pro- 
jected along each axis for / = / sat for a single run. This 
stronger measurement causes rapid decay of the condensate's 
spin. 



of the master equation over a single simulated experimen- 
tal run. We measure this magnetic field by observing the 
spin in a transverse direction; the atoms' Larmor pre- 
cession frequency reveals the field strength. We assume 
a condensate of 10 4 atoms confined within a transverse 
spatial extent of 15 /im in an ambient magnetic field of 1 
mG in the y-direction. We apply an additional magnetic 
field n — — in the negative z-direction; since G and f are 
well-controlled experimental parameters this should be 
possible to a high degree of accuracy. This additional 
field cancels out the classical portion of the light-matter 
interaction, though the photon shot noise fluctuations in 
Eq. (22) are still present. The probe is detuned by 27tx 
150 MHz below the F = 1 ->• F' = 2 (Dl) transition 
of 87 Rb, and interacts with the condensate for a total 
measurement time of 100 ms. The evolution of the con- 
densate spin is shown for probe intensities of (Fig. [lj 
I = 0.01/sat and (Fig. [2) I = J sat , where I sat is the satu- 
ration intensity. The measurement strength can best be 
characterized by the dimensionless ratio — i, where z/l 
is the Larmor frequency. For the parameters listed, this 
ratio is 0.1 and 10., respectively. In the former case, the 
free evolution of the atoms is not noticeably perturbed; 
yet the photocurrent signal unambiguously oscillates at 
the Larmor frequency. This will allow an accurate de- 
termination of the applied field. On the other hand, 
as the probe intensity is increased and the measurement 
strength exceeds unity, the backaction-induced stochastic 
evolution overwhelms the free Larmor precession, result- 
ing in a rapid decay of the transverse magnetization of 
the condensate in only a few oscillations. This crossover 
suggests a dynamical phase transition as one moves be- 
tween the weak and strong measurement regimes. 

Summary: We have provided a theoretical treatment 
of the quantum backaction due to the dispersive inter- 
action between a spinor Bose-Einstein condensate and 
an off-resonant light field. In addition to being the ba- 
sis for optical magnetometry using a Bose condensate, 
this interaction has also been shown to be a versatile 
quantum interface for quantum information processing 
and state engineering [15]. Straightforward additions to 
our model include a description of spatially inhomoge- 
neous spin textures in the condensate, stroboscopic opti- 
cal measurements of the condensate and the detection of 
time-varying magnetic fields. In addition to understand- 
ing the potential sensitivity of condensate-based mag- 
netic field sensing, our formalism can also be applied 
to quantum-limited nondestructive imaging of magnetic 
textures in spinor condensates and the creation of novel 
many-body states via quantum non-demolition (QND) 
measurement [16] . 
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